!  ad- fl 1 44  751 THE  EFFECTS  OF  RADIAL  YARNS  THREE -DIMENSIONALLY  1/  Q 

’  REINFORCED  CARBON-CARBON. . < U )  CALIFORNIA  UNIV  LOS 

ANGELES  SCHOOL  OF  ENGINEERING  AND  APPLIED.. 

UNCLASSIFIED  D  QUAN  ET  AL .  01  JUL  84  UCLA-ENG-84-22  F/G  11/4  NL 


MICROCOPY  RESOLUTION  TEST  CHART 


THE  EFFECTS  OF  RADIAL  YARNS  - 
THREEE  -DIMENSIONALLY  REINFORCED 
CARBON  -CARBON  COMPOSITES 

I ) V  i>m< 1  is  Qu.'in.  George  Sines,  S  3.  3  iMorf 


OIM R  MOD  14  ’7  C  050S 
I..  M  1  "  IliBS,  M  *n  :  ;r 


Unclassif led 


SECURITY  CLASSl*i«  ATloN  OF  THIS  PAGE  (Wlten  Date  htitered) 


REPORT  DOCUMENTATION  PAGE 


f’t  >H  T  NtjMHt  U 


r  I  I  l  t  >  an 1 1  Subtil  h< 

The  Kffects  of  Kadlal  Yarns  — 
Three -Dimensionally  Reinforced 
Carbon-Carbon  Composite 


& 


READiWSTRUCTIONS 
BEFORE  COMPLETING  FORM 


y  recipient's  catalog  number 


5  T  yPE  bF  REPORT  »  PERIOO  COVERED 

Technical  Report 


6  PERFORMING  ORG  REPORT  NUMBER 

UCLA- ENG-84 -2 2  / 


7  authors. 

Douglas  Quan,  George  Sines,  S.B.  Batdorf 

B  contract  or  grant  NUMBER^*) 

N0014-77-C-0505 

9  PERFORMING  ORGANIZATION  name  And  ADDRESS 

School  of  Engineering  and  Applied  Sciences 
University  of  California 

Los  Angeles,  CA  90024 

10.  program  element,  project,  task 

AREA  A  WORK  UNIT  NUMBERS 

H  controlling  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

12  REPORT  DATE 

1  July  84 

Att.  Dr.  A.M.  Dlness 

Arlington,  Virginia  22217 

13  NUMBER  OF  PAGES 

89 

14  MONITORING  AGENCY  NAME  ft  ADDRESSflf  differant  from  Controlling  Ollice) 

Monitor,  Dr.  L.H.  Peebles,  Jr. 

Office  of  Naval  Research 

Arlington,  Virginia  22217 

15.  SECURITY  CLASS,  (o t  thla  report) 

Unclassified 

15 a.  DECLASSIFICATION  DOWNGRADING 
SCHEOULE 

16  DISTRIBUTION  STATEMENT  (of  this  Report) 

Unlimited 

17  DISTRIBUTION  STATEMENT  (of  the  abatract  an  tar  ad  In  Block  20,  II  different  from  Report ) 

18  SUPPL EMENTARY  NOTES 

_ _ 

19  KEY  WORDS  ( Contlnu a  on  ravarae  aide  U  necaaamry  and  Identify  by  block  number) 


Carbon-Carbon  Composites 
Rocket  Nozzles 


Graphite 

Fabrication  Stresses 


70  ABSTRACT  (Continue  on  revera*  aide  If  necmaeary  and  Identify  by  block  number) 

Some  cylindrically  wound,  carbon-carbon  billets  have  had  gross  fracture 
of  the  circumferential  bundles  during  thermal  processing.  One  function  of 
the  radial  bundles  is  to  reduce  the  stress  that  causes  such  fractures.  An 
analysis  is  presented  to  show  the  potential  reduction  of  the  stress  in  the 
circumferential  bundles  during  processing  if  the  radial  bundles  remain 


DD  I  JAN *73  1473  EDITION  OF  1  MOV  63  IS  OBSOLETE 

S/N  0102  LF  014  6601  _ 

SECURITY  CLASSIFICATION  OF  THIS  FACE  Dm  It  tnltrmd) 


SECURITY  CLASSIFICATION  OF  this  PAae>'IF?l«n  Data  F.ntarad' 


intact.  A  simple  analysis  shows  that  the  stress  in  the  radial  bundles  is 
even  higher  than  that  in  the  circumf erent ials;  therefore,  they  are  likely 
to  fail  by  either  fracture  or  debonding.  The  radial  bundles  terminate  at 
the  outer  and  inner  radii;  an  analysis  is  made  of  the  debonding  from  the 
local  shear  near  the  ends  of  the  radial  bundles.  Partial  benefit  of  the 
radials  might  be  obtained  if  creep  reduced  the  stress  in  thorn.  An  experi¬ 
mental  study  on  the  creep  of  pitch- impregnated ,  uni-bundle  specimens  was 
conducted  and  results  are  presented.  Based  on  creep  of  the  radials,  a 
procedure  is  presented  to  find  an  optimum  time-temperature  path  to  avoid 
failure  of  the  radial  bundles  and  characteristic  optimum  paths  are  presented 
Exact  optimum  path  shape  cannot  yet  be  proposed  because  the  bond  strength 
and  bundle  strength  at  elevated  temperatures  as  well  as  the  creep  behavior 
are  not  yet  known  with  sufficient  accuracy. 


security  classification  of  this  PAOEfWi»n  lta-a  /  .  i 


The  Effects  of  Radial  Yarns, 


Three  Dimensionally  Reinforced  Carbon-Carbon  Composite 


Douglas  Quan,  George  Sines,  and  S.B.  Bntdorf 


Technical  Report 

for 

Contract  No.  0014-77-C-0S05 

Damage  Mechanisms  and  Modeling  of 
Carbon-Carbon  Composites 


July  1984 

Prepared  for 
Department  of  the  Navy 
Office  of  Naval  Research 
Arlington,  Virginia  22217 


i 

l 


>* 


A 


REPORT  ABSTRACT 


'Some  cylindrically  wound,  carbon-carbon  billets  have  had  gross  fracture 
of  the  circumferential  bundles  during  thermal  processing.  One  function  of 
the  radial  bundles  is  to  reduce  the  stress  that  causes  such  fractures.  An 
analysis  is  presented  to  show  the  potential  reduction  of  the  stress  in  the 
circumf erential  bundles  during  processing  if  the  radial  bundles  remain 
intact.  A  simple  analysis  shows  that  the  stress  in  the  radial  bundles  is 
even  higher  than  that  in  the  circumferentials; therefore, they  are  likely  to 
fail  by  either  fracture  or  debonding.  The  radial  bundles  terminate  at  the 
outer  and  inner  radii;  an  analysis  is  made  of  the  debonding  from  the  local 
shear  near  the  ends  of  the  radial  bundles.  Partial  benefit  of  the  radials 
might  be  obtained  if  creep  reduced  the  stress  in  them.  An  experimental 
study  on  the  creep  of  pitch- impregnated ,  uni-bundle  specimens  was  conducted 
and  results  are  presented.  Based  on  creep  of  the  radials,  a  procedure  is 
presented  to  find  an  optimum  time-temperature  path  to  avoid  failure  of  the 
radial  bundles  and  characteristic  optimum  paths  are  presented."*'  Exact 
optimum  path  shape  cannot  yet  be  proposed  because  the  bond  strength  and 
bundle  strength  at  elevated  temperatures  as  well  as  the  creep  behavior  are 
not  yet  known  with  sufficient  accuracy. 
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The  text  of  this  report  consists  of  the  masters  thesis  of  Douglas  Quan 
plus  additional  work  presented  in  four  addenda.  The  addenda  are  referenced 
in  the  body  of  the  report.  They  give  more  detailed  consideration  of  the 
uncertainty  in  the  measurement  in  this  study  of  the  linear  parameter  of 
the  creep  equation,  a  discussion  of  the  possible  physical  interpretation  of 
the  values  measured  for  the  activation  energy  for  creep  of  carbon-carbon 
composites,  analytical  demonstrations  of  the  effects  that  various  creep 
and  strength  parameters  would  have  on  the  optimum  time-temperature  path  to 
avoid  failure  of  the  radial  bundles,  and  some  more  detailed  consideration 


of  the  stress  in  the  radial  bundles. 


ABSTRACT  OF  THE  THESIS 


The  Effects  of  Radial  Yarns  in 
Three  Dimensionally  Reinforced  Carbon-Carbon  Composite 

by 

Douglas  C.  Quan 

Master  of  Science  in  Engineering 
University  of  California,  Los  Angeles,  1983 
Professor  George  H.  Sines,  Chairman 

Three-dimensional,  cyl indr ically  woven  carbon/carbon 
composite  is  particularly  appealing  for  rocket  nozzle  appli¬ 
cations  due  to  its  favorable  geometrical,  thermomechanical 
and  erosive  properties.  However,  the  performance  of  these 
materials  is  degraded  by  the  presence  of  various  kinds  of 
in-process  defects,  such  as  pores,  fiber  buckles  and  cracks. 
The  inhomogeneity  and  anisotropy  of  three-dimensionally 
woven  carbon/carbon  composite  prevents  accurate  modeling  of 
the  material  using  classical  elasticity;  therefore  elaborate 
computer  codes  for  the  finite  element  analysis  are  often  em¬ 
ployed.  Although  expensive  to  develop  and  employ,  these 
programs  are  an  indispensable  tool  for  the  analysis  of 
three-dimensionally  woven  composites.  However,  analytical 
treatments  using  simple  mathematics  can  provide  a  feeling 
for  the  physical  state  of  the  material  which  can  be  valuable 
to  the  designer.  The  intent  of  this  work  is  to  present  one 
such  analysis  using  a  simplified  approach. 
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M  =  mechanical 
T  =  thermal 
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[  nt  roduct  ion 


While  fracture  of  the  c i rcumf erent i al  yarn  during 
heating  in  fabrication  has  received  most  of  the  attention,  a 
simple  analysis  can  show  damage  to  the  radial  yarn  is  more 
likely  because  of  its  higher  stress.  The  well-being  of  the 
radial  fiber  is  very  important  to  the  integrity  of  the 
billet  as  a  whole.  When  the  billet  is  heated,  the  outer 
radius  will  be  forced  outward;  the  radial  yarn  wants  to 
prevent  this  from  happening  by  loading  itself  in  tension, 
thereby  relieving  some  of  the  load  in  the  hoop  yarns. 

Unlike  the  hoop  yarn,  the  radial  yarn  is  not 

continuous.  Both  ends  of  the  radial  yarn  intersect  a  free 
surface,  thus  a  zone  of  high  shear  must  be  established  at 
both  ends  of  the  yarn.  If  these  shear  stresses  are  greater 
than  the  yarn/matrix  interface  shear  strength,  the  fiber 
bundle  will  tear  away  from  the  surrounding  composite  so  that 
the  only  tensile  stress  in  the  bundle  are  those  resulting 
from  the  frictional  transfer  between  the  yarn  and  the 
surrounding  matrix.  Sudden  unloading,  such  as  debonding  of 
the  radial  yarn,  could  lead  to  catastrophic  failure  of  the 
circumferential  fiber  bundle  during  processing;  even  gradual 
unloading  can  lead  to  degradation  of  the  overall  strength  of 
the  finished  billet. 
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A  careful  crack  study  has  been  made  by  Kuhansedgh  and 
Sines  [1],  Some  of  the  findings  are  summarized  as  fellows: 

(1)  An  extensive  crack  structure  is  seen  in  the 

★ 

circumferential  surface  ,  which  in  some  cases  tends  to 

be  periodic;  that  is,  cracks  tend  to  occur  after  each 
five  or  six  circumferential  fiber  bundles. 

(2)  There  are  a  few  cracks  around  some  of  the  axial  fibers 
but  no  extensive  crack  pattern  is  present  around  these 
fibers  . 

(3)  Most  radial  fibers  were  totally  surrounded  by  well- 
defined  cracks.  Altho>ugh  some  appeared  to  have  we]  1  - 
defined  cracks  on  only  three  sides,  the  fourth  side  may 

also  have  been  debonded  and  these  cracks  were 
continuous  around  the  radial  fiber  tnroughout.  the 
billet ,  Fig .  1. 

Point  (3)  seems  to  support  our  claim  that  high  stress 

in  the  radial  fiber  bundle,  causes  the  debonding.  This  is 
because  if  the  crack  are  caused  only  by  shrinkage  the 
distribution  of  crack  should  be  quite  even,  since  the 

shrinking  of  the  matrix  material  is  assumed  to  be  isotropic. 
The  striking  contrast  in  the  finding  between  the  leas' 
stressed  yarn,  point  (2)  and  the  most  severely  stressed  yarn 

T  Ci rcumfer ential  cracks  are  those  in  the  surface 
perpendicular  to  the  radial  direction.  The  surface  contains 
the  circumferential  and  axial  directions. 


point  (3)  suggested  that  stress  level  in  the  fiber  must  play 
a  role  in  crack  formation. 

It  is  the  purpose  of  this  study  to  take  advantage  of 
the  creep  process  during  the  late  g raphi t i zat ion  cycles  and 
develop  a  heating  schedule  so  that  fiber  bundle  debonding 
will  be  minimized,  thereby  preserving  the  integrity  and 
minimizing  micro-cracking  which  is  not  caused  by  shrinkage 
o  f  the  matr i x  . 
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The  Billet  and  its  Fabrication 


The  NAVY-C4X1  billet  is  a  thr ee-d inons i o na 1  , 
cyl  ind  r  ical  1  y  woven  carbon/ca  rbon  composite.  The  bil1'*-  ?  ~ 
manufactured  by  General  Electric,  using  multiple  hint 
pressure  impregnation  process  at  15  ksi  with  OP277-15V  coal 
tar  pitch;  which  is  supplied  by  the  Allied  Chemicn' 
Corporation.  The  billet  is  woven  by  two  different  kinds  of 
yarns.  HM-10000  yarn  is  used  in  the  axial  and 
circumferential  directions,  and  the  HM-3000  yarn  in  the 
radial  direction.  10000  and  3000  corresponds  to  the  number 
of  filaments  in  each  yarn;  both  yarns  are  fabricated  by  the 
Hercules  Incorporated,  Magnamite  Graphite  Fibers  Division. 
The  variable  fiber  volume  fractions  are  shown  in  Fig."5,  and 
the  properties  of  the  fiber  and  pitch  materials  are 
presented  in  Table  (1)  and  Table  (2). 

Processing  of  the  billet  can  be  divided  into  two  steps: 
Automatic  weaving  and  densi f ication .  Fig. 3  shows  ;  he  process 
flow  chart  of  the  billet. 

Automatic  weaving  of  this  thr ee-d imens i onal 1 v 
reinforced,  carbon/carbon  composite  billet  consisted  of 
fabricating  and  assembly  of  yarn  bundles  of  three  orthogonal 
directions  (radial,  axial,  circumferential).  The  billot  is 
then  heated  to  a  minimum  temperature  »f  150  C  (552  F )  before 
the  densi f ication  process  begins.  In  the  dousi f i ca t ion 
process  the  billet  was  densi fied  by  subsequent  impregnation. 


a 


carbonization,  and  g r aph i t i za t  io n  steps  using  pitch 
materials  as  the  impregnant.  During  the  carbonization  stage 
of  the  first  cycle  the  matrix  material  became  g raph i t i zab le 
carbon,  forming  mesophase  during  the  heating  treatment  in 
the  range  from  400  C  (750  F)  to  450  C  (840  F) .  Mesophase 
consists  of  large  po 1 ycond ensed  aromatic  compounds  oriented 
nearly  parallel  to  each  other  but  still  remaining  in  the 
] iquid  state  . 

Impregnation  of  carbon/carbon  composite  materials 
started  after  the  first  cycle  of  carbonization  and 
graphi tization  in  order  to  fill  up  the  cracks  which  were 

generated  during  those  steps.  This  step  known  as  vacuum 
impregnation,  was  carried  out  at  240  C  (480  F)  under  a 
pressure  of  40  torr.  The  degree  of  impregnation  depends  upon 
the  processing  conditions  such  as  pore  size,  pore  size 
distribution,  and  gas  permeability. 

High  pressure  carbonization  was  the  second  stage  of  the 
densi f ication  process.  The  billet  was  heated  to  550  C  (1200 
F)  where  the  material  carbonized;  some  cracking  due  to  this 
phase  change  resulted  in  a  relatively  stress-free  material 
at  this  step.  The  amount  of  cracking  at  this  stage  is 


observed  to  be  small  compared  to  cracking  in  subsequent 
steps  . 


The  last  step  is  g r aph i t i za t io n ,  the  billet 


wi  s 


subjected  to  a  uniform  temperature  rise  up  to  27sn  g  (V'nn 
F) .  At  this  stage  some  chemical  and  structural  changes 
occurred  in  the  matrix.  The  densi f ication  process  was 
repeated  four  times,  the  latter  three  in  order  to  replace 
material  volatilized  during  carbonization  and  g r aph i t i za t io n 
as  well  as  to  fill  the  cracks  occurring  each  time  on  cooling 
from  graph i t i za t io n .  The  billet  was  then  subjected  to 
inspection  and  machining. 


l_.  Estimation  o f  Stress  i n  t  ho  Radial  Bundle 

The  following  analysis  gives  an  estimation  of  the 
relative  magnitude  of  the  stress  in  the  radial  bundle  to 
that  of  the  circumferential  bundle.  For  the  purpose  of 
simplicity  we  assumed  that  the  circumferential  stress  at  the 
outer  and  inner  radius  are  equal  and  opposite. 
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which  implies  that  (  )Q  =  -(  =  I  1  From 
studying  the  results  obtained  by  [1]  and  [3],  this  found  to 
be  an  acceptable  approximation.  We  also  assumed  that  the 
Young's  modulus  in  the  direction  along  the  fiber's  axis  are 
equal  for  both  radial  and  circumferential  bundles  so  that 

(  Err  )a  =  (  )a  =  E.  Table  (1)  show  that  the  difference 
in  moduli  between  these  two  types  of  fiber  bundles  in  the 
navy  billet  is  about  eight  percent.  Lastly  the  effect  of  the 
axial  bundles  were  not  considered  in  this  analysis. 

The  circumference  of  the  billet  is  C  =  2"ir"R  When  the 
billet  is  heated,  the  inner  and  outer  circumferences  of  the 
billet  will  be  forced  to  change  its  dimensions.  The  outer 
circumference  just  after  a  rapid  heat  up  can  be  expressed 
as : 

I 

C  =  C  +  A  C  =  C  +  C  <  +  C  oAT 

O  0^0  0  O 


where  <  is  the  mechanical  strain. 


2"Ro  =  2FRd  +  2FR0!-|^j  +  2¥RooAT 

d-n' 

»:  -  *. [» ♦  ^  ^  i 

( I  -  3 ) 

where  oc  is  the  coefficient  of  thermal  expansion,  R  is  the 

f 

original  billet  radius  and  R  is  the  billet  radius  just 
after  a  rapid  heat  up. 


Similarly,  the  inner  radius  will  be. 


To-  i 

Ri  -  Ri  1  -  —  +  oe^r  j 


'’’he  mechanical  strain  in  the  radial  direction  is 


(1-4) 


rr 


[<Ro  -  V  -  ">0 

1 

-  VI 

1  "o  -  Ri 

J 

ocAT 


(I-r  ) 


Tr  +  R.l  ov 
.  _  j_o _ i_{  Re 

rr  i  R  -  R .  |  e 


(I-h) 


This  is  the  maximum  mechanical  strain  that  a  radial  bundle 
will  experience,  if  perfect  bonding  exists  between  the 
radial  and  the  boundaries  circumferential  bundles.  Stress  in 
the  radial  bundle  would  be: 


cr  =  E  < 
r  r  r  r  r  r 


rr  IE 


ee 


+  R.l 

O 

1  1 

R 

-R.l 

o 

'  J 

(  T  -■») 
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Since  we  have  assumed  E 


along  the  fiber's  axis 


rr  '  Eee 


^rr  =  *90 


(I-R* 


From  the  above  analysis  one  can  see  that  the  stress  in  the 

radial  bundles  is  always  greater  than  that  of  the 

circumferential  bundles.  For  a  billet  with  short  radial 

fiber  or  small  R  -r.  relative  to  the  billet  circumference 

o  i 

will  result  in  very  high  radial  stress,  which  agrees  with 
physical  intuition. 


*  Addendum  IV  further  discusses  the  approximations  made  in 
this  analysis.  It  also  gives  the  results  of  a  second- 
order  approximate  analysis  made  by  Julius  Jortner. 
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I_-A  The  Restraining  Effect  of  the  Rail  i  al  Bundle 

The  effect  of  the  radial  bundle  on  the  billet's 
circumferential  bundles  can  be  demonstrated  by  modifying  the 
computer  model  developed  by  Kuhansedgh  and  Sines  fll.  Tn  ( 1 1 
the  billet  was  modeled  as  an  assembly  of  concentric  thin 
walled  cylinders.  Alternate  cylinders  have  different 
properties;  one  set  models  the  fiber,  the  other  the  matrix. 
The  number  of  cylinders  modeling  the  fibers  are  equal  to  the 
number  of  layers  of  circumferential  wrap.  The  Navy  billet 
used  in  this  study  has  59  cylinders  (  numbered  from  1  r0  5B 
).  The  odd  numbers  are  for  the  fiber  cylinders  and  the  even 
numbers  the  matrix  cylinders.  Each  cylinder  can  be  expressed 
in  a  general  form. 


( I A  —  1  *) 


where  n=l,2,3,...  59 

When  n  =  odd  E  =  modulus  of  fiber 
x 

E  =  modulus  of  matrix 

y 

n  =  even  E  =  modulus  o f  matrix 
x 

Ey  =  modulus  of  fiber 

The  model  also  assumed  that,  the  properties  ,,i  •  h,>  fiber 


and  the  matrix  cylinders  will  not  change  with  temperature 


.md  Mu1  cono-ntr  ic  cylinders  are  ini  tally  in  stress-free 
contact . 


A  change  in  temperature  causes  the  cylinders'  diameters 

to  change,  and  thereby  introduce  i n ter -cyl ind r i cal  pressure 

P,  ...  throughout  the  billet  radii.  These  inter¬ 
in, n+1) 

cylindrical  pressures  can  be  obtained  by  solving  a 
simultaneous  set  of  equations  which  takes  the  form  of 
equation  (  IA-1  ) ,  these  pressures  can  then  be  used  to 
calculate  the  value  of  the  billet's  hoop  fiber  stress  with 
the  radial  effect  omitted  Fig. A. 


The  restraining  effect  of  the  radial  fiber  can  be 
modeled  by  applying  a  hydrostatic  pressure  at  the  inner  and 
outer  surfaces  of  the  cylindrical  walls.  The  magnitude  of 
the  hydrostatic  pressure  is  equal  to  the  traction  induced  by 
the  highly  stressed  radial  fiber  on  a  characteristic  area 
Fig. 5.  Stress  on  the  radial  fiber  and  its  cross  sectional 
area  is  assumed  to  be  constant  along  its  length,  but  the 

characteristic  area  decreases  as  you  move  towards  the  center 
of  the  billet;  one  should  immediately  realize  that  the 
hydrostatic  pressure  acting  on  the  inner  surface  is  greater 
than  the  pressure  at  the  outer  surface.  Note  that  these 
hydrostatic  pressures  represent  the  restraining  effect  of 
the  radial  bundles  on  the  hoop  bundles,  therefore  we  can 
only  consider  the  amount  of  stress  that  has  successfully 
transferred  to  the  hoop  bundle  in  calculating  the 
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hydrostatic  pressures.  In  this  work  the  load  transfer 
efficiency  was  assumed  to  be  0.50,  it  is  a  reasonable 
assumption  considering  the  amount  of  debonding  and  cracks, 
which  is  so  commonly  fc'und  around  a  fiber  which  intersects  a 

free  surface. 

To  calculate  these  restraining  pressures  the  original 
billet  model  (with  no  radial  effect)  was  Jtilized  to  find 
the  stress  of  its  hoop  fibers.  Knowing  these  stresses  one 
can  then  calculate  the  change  in  billet  radius  just  after  a 
rapid  temperature  increase.  Now  a  fictitious  radial  fiber  is. 
super  imposed  onto  the  billet  with  its  ends  fixed  to  rh<' 
outer  and  inner  most  hoop  fibers,  one  can  then  estimate  the 
radial  stress  induced  by  the  radius  change.  RecaV,  an 
assumption  has  been  made  that  only  50%  of  the  stress  was 
successfully  transferred  and  working  to  keep  the  billet  fren 
expanding,  therefore  the  force  in  the  radial  bundle  is 
multiply  by  a  factor  of  0.5.  divide  these  force  by  their 
characteristic  area  and  the  results  are  the  hydrostatic 
pressures.  See  appendix  for  detail  mathematical 
illustrations . 

The  effect  of  these  pressures  can  be  seen  by 
substituting  them  into  the  billet  model.  The  results  are 
plotted  on  Fig. 6.  This  plot  has  over  estimated  the 

effectiveness  of  the  radial  bundle  from  a  real  bil  let  .  Th  i s 
is  because  the  radial  stress  which  was  used  h  calculate  the 
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hydrostatic  pressures  was  computed  from  a  billet  where  the 
maximum  possible  amount  of  expansion  was  allowed.  In  actual 
case  internal  frictional  forces,  anisotropy  stresses  and 

stuffing  stresses  will  prevent  the  hoop  fiber  from  expanding 
this  much.  Nevertheless  the  usefulness  of  a  strong  and  well 

bonded  radial  fiber  is  reviewed.  In  this  model,  the 

presence  of  the  radial  fibers  eliminated  the  entire 
compressive  stress  region  and  decreased  the  stress  level  of 
the  outer  portion  of  the  billet  by  more  then  20%.  The  sharp 
drops  in  stress  level  that  you  see  in  this  plot  is  due  to 
the  model's  extreme  sensitivity  to  the  variation  of  fiber 
volume  fraction  Fig. 2. 

In  an  actual  billet  a  strong  and  well-bonded  radial 
yarn  is  expected  to  help  to  reduce  some  of  the  processing 
anomalies.  A  smaller  compressive  zone  can  reduce  mini- 
buckling  at  the  inside  diameter,  lower  tension  at  the  outer 

diameter  can  reduce  the  chance  of  cracks  formation  in  the  rz 
plane.  One  would  expect  a  flatter  and  smoother  hoop  stress 

distribution  for  the  actual  billet.  In  this  analysis  it  was 
found  that  58%  of  the  hoop  fibers  have  a  stress  value  1/3  of 
that  in  the  radial  fibers  and  42%  of  the  the  hoop  fibers 
have  a  stress  value  1/2  of  that  in  the  radial  fibers.  These 
figures  again  reinforced  our  claim  that  the  radial  yarn  is 
the  most  highly  stressed  bundles  in  the  billet  and  its 

well-being  is  crucial  to  the  integrity  of  the  billet  as  a 
whole  . 
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1 1 .  Fiber/Matr ix  Displacement  Mismatch  Analysis 


Property  mismatch  between  the  fiber  and  matrix  caus : ng 
thermal  stresses  could  lead  to  the  ultimate  failure  of  the 
cyl indr ically  woven  carbon/carbon  composites  durino 
fabrication.  Adhesive  failure  at  the  radial  fiber  cr  failure 
of  the  fiber  itself  is  expected  to  take  place  before  any 
other  mode  of  failure.  This  is  because  the  radial  yarns  are 

the  most  highly  stressed,  and  regions  of  high  shear  stress 
are  created  at  the  free  ends  because  of  property  mismatch. 

Imagine  a  billet  composed  of  equally  spaced  fibers  with 
matrix  in  between,  as  shown  in  Fig. 7.  When  heated,  the 
fibers  and  matrix  will  expand,  but  the  matrix  will  try  to 
expand  more  because  of  its  higher  coefficient  of  expansion. 

This  mismatch  will  result  in  high  shear  stress  at  the 
fibers'  ends.  If  the  stress  level  exceeded  the  bond 
strength  of  the  interface,  debonding  will  initiate.  Using 
some  physical  intuition,  one  can  realize  that  the  debcnded 

portion  of  the  matrix  will  try  to  flatten  itself  out,  which 

translates  into  creating  a  pulling  force  on  the  matrix 
around  the  matrix-fiber  interface  Fig. 8.  If  this  situation 

were  to  occur,  the  usefulness  of  the  radial  yarn  in  terms  of 
load  carrying  capacity  is  reduced  by  the  debonding. 

The  fiber-matrix  displacement  mismatch  can  be  analyzed, 
utilizing  a  model  very  similar  to  the  one  discussed  earlier 
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in  this  section,  Fig.  7. 

Balance  of  force  exist  between  the  fiber  and  the  matrix 
during  the  heating  up  process  is  expressed  by, 


M  M 

Af  Ef  =  Am  E 
t  t  t  m  m  m 


(II-l) 


The  difference  in  mechanical  strain  is  induced  by  the 
mismatch  in  coefficient  of  expansion: 


M 

f 


(oc  -  ocF)  (T  -  T1) 


(II-?) 


where  is  the  initial  temperature 

Let  S  be  the  shear  force  between  the  interacting  surfaces, 
and  has  a  dimension  of  force  per  unit  length. 


S  =  kj  (Uf 


U  ) 
m ' 


where  is  a  constant  and  U  is  displacement. 


(11-3' 


The  same  shear  force  can  be  written  in  terms  of  displacement 
of  the  fiber,  or  displacement  of  the  matrix,  which  results 
from  their  respective  strain. 


S 


d  U  r 

d(AfRfir>r) 

3x 


If  and  E^  are  constant, 


fII-4) 
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differential  equation  of  the  form 


d. -r?  -  k?d  =  n 


(ir-io) 


2  2 

The  character istic  equation  is  X  -K‘=0,  which  can  be 
factored  into  (K~K)(X+K)  =  0.  The  roots  are  X^=K 
Since  the  roots  are  real  and  distinct,  the  solution  is 


n  Kx  -Kx 

D  =  c^e  +  Cpe 


(11-11) 


For  K=-K  we  can  write 


D  =  kjCosh  Kx  +  k^sinh  Kx 


(11-12) 


kl  =  cl+c2  k  2  =  cl~c2 


(11-13) 


In  obtaining  the  above  solution,  we  use 


cosh  Kx  =  Y(®KX+e  KX) 


<II-H) 


sinh  Kx  =  i(eKx-e  Kx) 


The  boundary  conditions  are: 


UfO)  =  0 


(I  I -1 S ) 


U’  (|) 


(11-16) 


for  U(0)=0 
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0=ci+Cp  cl=-c2 


for  U*  (^)=  (<xm-ocf  )^T 


(am-af)^T  =  CjKe^-c^Ke  Kx 


(II-1 


=  c1KIeKx+e  Kx) 


=  2c^K  cosh  Kx 


(«  "«f)AT 

c^  =  - -  cosh  Kx 


( 1 1  - 1  Q  ' 


D  =  ci(eKx-e  Kx) 


=  c^2  sinh  Kx 


(«m-<xf)AT  sinh  Kx  («m-«£)Ar 


K  cosh  Kx 


Tanh  Kx 


( T I -1 9  > 


Note  that  the  total  displacement  mismatch,  n  is  zero  at 
the  center  and  is  maximum  at  the  ends,  equation  (TT— 19) . 
Large  mismatch  means  great  shear  stress  at  the  ends,  which 
encourages  the  initiation  of  adhesive  failure. 


III.  Critical  Tensile  Stress 


In  the  previous  section,  it  was  shown  that  the  property 
mismatch  between  the  fiber  and  the  matrix  creates  a 
displacement  mismatch.  For  the  following  analysis 
everything  surrounding  a  radial  fiber  will  be  considered 
matrix.  This  assumption  can  be  argued  by  looking  at  an  unit 
cell  of  the  material.  Fig. 9.  The  elastic  modulus  of  the 
axial  and  circumferential  fibers  are  low  along  the  radial 
direction.  At  elevated  temperature,  the  transverse  modulus 
of  the  fibers  are  very  close  to  that  of  the  matrix, 
reference  [4].  In  the  light  of  this  approximation,  the 
assumption  that  the  radial  fiber  is  completely  surrounded  by 
matrix  for  stress  analysis  of  the  radial  fiber  along  the 
radial  direction  will  be  valid. 

The  shear  stress  induced  by  the  displacement  mismatch 
caused  by  the  difference  in  thermal  expansion  coefficient  is 
analogous  to  the  shear  stress  induced  by  displacement 
mismatch  caused  by  the  difference  in  elastic  modulus.  The 

latter  problem  has  been  treated  by  Cox,  [51,  using  the  so 
called  "shear  lag  analysis".  In  this  model  a  mechanical 
stress  is  applied  on  the  infinite  resin  parallel  to  the 
embedded  yarn.  In  the  region  of  the  yarn  ends  the  strain  of 
the  fiber  will  be  less  than  the  average  strain  in  the 
matrix.  The  deformation  field  in  the  matrix  has  been 
sketched  in  Fig. 10.  Cox  shown  that  the  tensile  stress  along 


this  yarn  is  given  by: 


ov  = 


Ef 

t  m 


,  i 

1-cosh  p  f^d-x)  I 


cosh  -^3d 


4 


(  T  T  T -1  ) 


where 


P  - 


2G 


m 


Efr2ln(i) 


2L 

is 

the 

inter  f iber 

spac ing 

r 

is 

the 

radius  of 

the 

fiber 

d 

is 

the 

length  of 

the 

fiber 

and  the  shear  stress  at  the  interface  is  give  by: 

< 


r  = 


Ef«m 
t  m 


m 


PEf  1  n  ( 


|2sinh  pf^d-x) 


1 

IS' 


n  it 


f  T  I  T  -  M 


The  ratio  of  maximum  interface  shear  stress  at  the  fiber 
ends  to  maximum  tensile  stress  in  the  fiber  can  be 
calculated  by  equations  ( 1 1  — 1 >  -  (II-3),  and  is  found  to  he: 


max 


7» 


f '  max 


coth  ii™ 


2Rf  Inf—) 


flTI-A) 
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r  f  (■  he  max  iinum  shear  si  r**:::;  that  in  interface  ran  withstand 
can  be  estimated,  either  experimentally  or  analytically, 
then  we  can  calculate  the  tensile  stress  in  the  fiber  bundle 
using  equation  (III— 4)  .  Knowing  this  critical  stress  value, 
cr*  we  can  then  tailor  our  heating  schedule  so  that  the 
stress  in  the  radial  fiber  will  never  exceed  this  critical 
value,  thereby  minimizing  micro-cracks  that  are  caused  by 
adhesive  failure  due  to  displacement  mismatch  during  the 
heating  up  process. 

It  should  be  noted  that  the  analysis  used  here  is  not 
exact.  Finite  element  analysis  and  experimental  studies 
suggested  that  it  underestimates  the  shear  stress 
concentration  at  the  ends  of  the  fiber  by  about  a  factor 
two.  One  should  always  adjust  the  results  from  equation 
( 1 1 1  — 3 )  by  a  factor  of  at  least  two  for  application,  [51. 
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IV.  Heating  Rate  -  Stra in  Mismatch  Approach 

During  the  manufactur i ng  of  the  billet,  the  material  is 
subjected  to  a  series  of  heating  and  cooling  cycles.  The 
significance  of  shrinkage  cracks  caused  by  the  contraction 
of  matrix  material  during  cooling  is  undeniable.  Many  oeople 
have  overlooked  the  possibility  that  micro-cracks  can  be 
generated  by  displacement  mismatch  between  the  fibr>r  ^ n<> 
matrix  during  heating  up  for  reasons  that  has  been  discussed 
in  previous  sections. 

Recall  it  has  been  mentioned  in  section  IT  that  the 
mechanical  strain  mismatch  can  be  expressed  in  terms  •, f 
thermal  expansion  coefficient  mismatch; 


=  (ocm-ocf)dT 

Then  the  strain  rate  is 


d< 

cTF 


,  ,dT 

^m-^cft 


( I  v  - 1  ) 


(TV-?  ) 


The  assumption  has  been  made  that  the  billet  was  subjected 
to  a  rapid  heating  process,  so  that  no  significant  creen  can 
take  place  during  the  heat  up.  The  tensile  stress  in  the 
yarn  at  this  point  will  therefore  independent  of  time.  By 
the  same  token  the  balancing  compressive  stress  in  the 
matrix  must  also  be  time  independent.  An  excessively  quick 
heating  rate  is  undesirable,  since  it  would  lead  to  adhesive 
failure  at  the  fiber  matrix  interface.  A  controlled  he.nt:ni 
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rate  must  be  used  in  order  to  allow  creep  to  partially 
balance  out  the  strain  caused  by  the  mismatch.  Assuming 
creep  only  take  place  in  the  yarn,  the  strain  mismatch  which 
is  to  be  relieved  by  the  creep  of  the  fiber  can  be  equated 

to  the  general  creep  equation; 


d<f  _ 

<TF  " 


^  ( o* )  n 


expf-^ipl 


( I V-3 ) 


it  , 

cr  =  is  the  critical  tensile  stress  in  which  a  fiber 
can  withstand  without  debonding  to  occur 
A  =  material  constant 
n  =  stress  exponent 
=  activation  energy 

R  =  gas  constant 
T  =  temperature 


If  oc,  E,  Aft  r  a  and  n  are  all  temperature  independent, 


dT 
d  t 


r^ii 

RT 


fIV-4) 


oc  -ocf 
m  t 


A(cr  ) 


n  J  exP  dT 


CIV-5 ) 


If  all  the  parameters  are  known,  then  the  rate  of 
temperature  change  can  be  found  with  the  tensile  stress  in 


the  fiber  holding  at  a  desirable  constant  value.  The 
integral  of  equation  (IV-5)  will  have  to  be  evaluated 
numerically.  The  resulting  plots  will  take  the  shape  shown 
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schematically  in  Fig. 11a,  with  the  stress  vs.  time  cur>/e 
holding  constant  as  shown  in  Fig. lib. 
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V.  Heating  Rate  -  S tress  Relaxation  Approach 

When  a  billet  is  heated  to  a  elevated  temperature,  one 
would  expect  the  highly  stressed  radial  fiber  will  try  to 
stretch  more  along  the  radial  direction  than  the  hoop  fiber 
would  allow.  If  a  strong  bond  exists  between  the  radial  and 
circumferential  fiber,  this  will  be  analogous  to  the  classic 
model  of  stress  relaxation  where  the  stress  in  a  heated 
fixed  ends  rod  is  being  relaxed  due  to  creep. 


Assume  that  the  only  meaningful  creep  action  is 
performed  by  the  radial  yarn.  Recall  the  unit  cell  of  Fig.R, 
where  the  assumption  was  made  that  25%  of  the  cell  is  radial 
fiber  and  the  remaining  75%  has  properties  similar  to  matrix 
along  the  radial  direction.  If  balance  of  force  exit 
between  the  matrix  and  fiber. 


* 

A 

m 


cr 

m 


AfO-f 


where  A  =  effective  area  of  matrix, 
m 


0. 75  o-  =  0. 25  cr 
m 


(V-l) 


(V-2) 


rv-2 ) 


Since  the  stress  term  in  the  creep  equation  behaves  in  a 

power  law  fashion,  this  makes  the  matrix  in  the  radial 
direction  relatively  insensitive  to  creep. 
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For  the  creep  of  the  redial  fiber,  use  the 


creep  equation: 


d< 

ar  ~ 


Ac- 


ex  p 


RT 


where  <  =  creep  strain. 


Since  the  creep  curve  of  carbon  fibers  are  linear,  therefor-' 
it  is  appropriate  to  apply  the  steady  creep  approximation  - 
elastic  analogy. 


ldo  n 

E  dt  '  A<r 

ex  p 

-AH  _  n 

TFr  -  n 

integrating  both 

sides  g ives 

1-n 

Z -  =  AF.t 

1-n 

exp 

lAH  +  c 

RT 

(V-h ' 

Use  the  initial 

condition 

to 

evaluate  the 

i n teg  ration 

constant 


At  t=0 ,  cr  =  <>e 

where  o-  =  stress  cause  by  elastic  strain  during  rapid  heat 
up,  before  any  significant  creep  has  taken  place. 

c  is  found  to  be 


c 


1  -n 


1  -n 


(V-"7  1 


Substituting  back  into  eqn  { V— f> ) 


1  -n 


cr  =  E-* 
e  e 


where  <tq  =  elastic  strain 


(V-10) 


cr  =  £(l-n)AEt  exp  +  (Ecrj  I_nj^  " 


CV-1 1 ) 


A,  n,  E,  Aji  are  to  be  determined  by  experiment.  For  carbon 
fibers,  n  has  a  value  greater  than  one.  Equation  (V-1.1) 
describe  the  stress  in  a  fiber  at  any  time  t  for  a  given 
temperature  T. 


The  purpose  of  regulating  the  stress  in  the  radial 

fiber  is  to  avoid  having  the  stress  exceed  a  critical  stress 

* 

value  <5  .  The  temperature  versus  time  and  the  corresponding 
stress  versus  time  plot  will  have  shape  such  as  those  shown 
in  Fig. 12.  This  heating  schedule  may  be  more  practical  in 
actual  manufacturing  practice;  unlike  the  previous  method  in 
which  the  temperature  controller  has  to  follow  the  exact 
path  described  by  equation  (IV-5),  this  approach  allows  an 
approximation  to  that  path. 
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The  stress  during  the  rapid  heat  up  < eo  .  a-h  of 
Fig. 12b)  can  be  calculated  quite  easily,  because  the 
assumption  has  been  made  that  durinq  the  heating  process  no 
creep  will  occur.  Therefore  stress  is  caused  by  both  elasti~ 
and  thermal  strain.  When  the  stress  value  reaches  point  b, 
the  critical  stress  value,  temperature  should  hold  constant 
(  such  as  b'-c’at  Fig. 12a  )  so  that  stress  can  be  relayed 
(points  b-c  in  Fig. 12b)  before  further  heatinq.  The  holdioi 
time  has  to  be  decided  by  the  user  depends  on  the  t yoe  of 
fiber.  The  holding  time  should  not  be  too  long  for  obvious 
economic  reasons. 


V I .  The  Creep  Rehav i o r  o f  an  Impregnated  Graphite  Fiber 
Bundle 


The  objective  of  this  study  is  to  determine  the  stress 
dependence  n  and  the  activation  energy  H  terms  of  the  well 
known  empirical  creep  equation. 


d  4 
TtT  “ 


Act 


ex  p 


rVT-1) 


of  an  impregnated  graphite  fiber  bundle. 

The  graphite  fiber  used  in  this  study  is  Magnamite 
graphite  fiber  type  HM-3000  which  is  a  continuous  high 
modulus,  PAN-based  fiber  with  3000  filament-count  tows.  The 
impregnant  used  is  the  CP277-15V  coal  tar  pitch  manufactured 
by  the  Allied  Chemical  Corporation. 


VI -A  Specimen  Preparation 

Before  the  specimen  can  be  tested  in  the  high 
temperature  creep  furnace,  a  series  of  task  will  have  to  go 
into  preparing  the  specimens.  A  half  inch  diameter  loop  was 
formed  at  both  ends  of  the  dry  yarn  bundle  and  was  secured 
with  twelve  knobs  distributed  over  a  two  inch  section. 
Fig.  13.  Cotton  sewing  tread  was  used  for  making  the  knobs, 
they  should  be  strong  enough  to  prevent  any  slipping  while 
being  pulled  at  both  ends.  The  dry  yarn  measures  16  inches 
from  end  to  end  and  is  now  ready  for  impregnation. 
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The  15V  coal  tar  pitch  was  pretreatod  at  300C  under  a 
nitrogen  atmosphere  for  about  five  hours  before  it  was  used 
in  the  actual  impregnation  of  the  fiber.  The  reason  fc r  this 
pretreat  process  was  to  drive  out  some  of  the  undesirable 
chemicals  such  as  sulfur,  ash  and  hydrogen  which  made  up  the 
pitch.  A  typical  coal  tar  pitch  has  a  carbon  content  of  93?:, 
a  hydrogen  content  of  5.5%,  a  sulfur  content  of  1.45%  and  an 
ash  content  of  0.05%.  These  chemicals  are  highly  volatile  at 

high  temperature  therefore,  using  untreated  pitch  as 
impregnant  can  cause  problems  such  as  the  building  up  c.  f 
sulfur  and  ash  particles  causing  the  blockage  of  the 
nitrogen  gas  exit  tubing  and  the  bloating  of  the  impregnated 
yarn  during  the  later  calcination  process. 

The  graphite  fiber  bundles  were  impregnated  in 
pairs, they  were  laid  into  a  18  x  .75  x  .5  inch  trench  cut 

out  from  a  rectangular  cross-sectioned  aluminum.  A  highly 
sanded  wooden  stick  was  placed  at  each  ends  through  the 
hoops  of  the  fiber  bundles  Fig.  14.  The  function  of  the 
sticks  was  to  keep  the  ends  of  the  fibers  above  the  pitch 
and  for  stretching  the  bundles  when  it  was  being  raised  from 
the  resin  bath.  The  trench  was  then  filled  with  pretreated 

15V  pitch  powder  and  was  heated  on  a  hot  plate.  The  pitch 
melts  at  around  HOC  and  the  bundles  were  allowed  to  sit  in 
the  resin  bath  for  about  2.5  hours  to  ensure  good  wetting 
throughout . 
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300C 
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wa  s 

100C/hr.;  :from  300C  to  750C  which  is  the  critical 
temperature  range  the  heating  rate  was  about  15C/hr.;  from 

75 OC  to  900C  the  heating  rate  was  about  50C/hr.  The  heating 
history  of  each  pair  of  fiber  bundles  was  recorded  on  a 
plotter  for  later  reference.  At  about  900C  the  furnace  was 
turned  off  and  allow  to  cool  to  room  temperature  with 
nitrogen  flowing  at  all  time  to  insure  minimal  oxidation. 
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VI -B  Creep  Exper iment ;  Instrumentations  and  methodology 

The  creep  experiment  was  performed  at  a  facility  in  the 
Aerospace  Corporation's  Material  Science  Laboratory.  Tht 
instrumentations  used  are  shown  in  Fig. 19  ana  2D.  The 
calcined  fiber  bundle  was  held  in  place  inside  the  hiuh 
temperature  creep  furnace  by  two  mechanisms  shown  in  Fig. 2’, 
These  clamping  mechanisms  were  machined  out  from  a  very  Cine 

grained  graphite.  One  end  of  the  fiber  was  fixed  and  tie 
other  end  unrestrained,  the  unrestrained  end  was  connected 

to  a  weight  which  passed  over  a  smooth  pulley.  A  linear 
variable  digital  transducer  (LVDT )  is  also  connected  to  this 
free  end,  any  yarn  displacement  would  be  detected  by  the 
LVDT.  The  signal  sent  out  by  the  LVDT  will  then  be  measured 
by  a  digital  volt  meter  (DVM)  and  recorded  in  voltage  by  a 
HP9815  computer.  The  computer  was  so  programmed  that  a 
reading  was  taken  every  three  minutes  and  forty-five 
seconds.  When  the  computer  takes  a  reading,  a  glitch  maker 
which  was  connected  to  the  computer  and  a  two  channel  chart 
recorder  would  make  a  spike  on  the  chart  at  tha*-  instant 

which  aids  in  later  data  interpretations.  The  furnace 
current  was  controlled  manually  below  1500C  and  was  then 

switched  over  to  computer  monitoring  over  1S00C.  Below 

1500C  the  temperature  was  measured  by  a  manual  optical 
pyrometer,  at  above  1500C  the  temperature  was  measured  by  a 

IRC  ON  Infrared  pyrometer  which  was  connected  to  the  sane  twc 
channel  chart  recorder,  which  means  that  the  chart  recorder 
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will  not  start  recording  temperature  readings  until  the 
furnaced  reached  1SOOC.  Since  significant  creep  usually  take 

Pi  ace  at  temperature  above  2200C  for  this  material,  no 
important  data  is  lost  from  such  an  arrangement. 


The  stress  dependence  of  the  creep  rate  was  determined 
by  a  change-of-stress  method  in  which  the  structural  changes 
which  usually  take  place  during  creep  is  minimized.  In  this 
unispecimen  method  each  specimen  was  strained  at  a  selected 
temperature  under  an  argon  atmosphere  and  stressed  until  the 

creep  curve  was  well  established,  at  which  time  the  stress 

was  changed  suddenly  to  another  value  while  holding  the 
temperature  constant.  The  ratio  of  the  creep  strain  rate 

immediately  after  the  stress  change  to  that  just  before  is 
compared  to  the  ratio  of  stress  after  the  change  to  that 
immediately  before  to  determine  the  stress  dependence  of  the 
creep  rate.  The  stress  dependence  was  found  to  take  the 
following  form: 


n 


•  • 

In  < 2// * 
In  ct^/o- 


1 

1 


(VI-2) 


where  and  are  the  creep  rates  immediately  after  and 
immediately  before  the  change,  and  Co,  and  are  the 

respective  stresses.  Since  the  structure  immediately  before 
and  after  the  stress  change  was  essentially  the  same,  the 

structural  differences  that  normally  complicates  analyses  of 
the  stress  dependence  of  the  creep  rate  does  not  exists. 
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The  same  line  of  reasoning  is  applied  to  determine  the 


activation  energy  of  the  impregnated  bundle.  This  tim<>  we 
change  the  temperature  while  holding  other  variables,  stress 
and  structure  constant.  Knowing  the  temperature  difference 
and  the  creep  rate  at  respective  temperatures,  the 
activation  energy  is  found  to  be  in  the  following  form: 


R  Ti  T2r 

^  *  t2  -  Tx[  ln  42 


In 


fVT 


The  scattering  of  and  n  values  is  expected  despite 
the  outmost  care  that  has  gone  into  preparing  and  testing 
the  specimens.  A  best  fitted  value  will  have  to  be  chosen 
for  each  of  the  two  parameters,  they  are  in  turn  used  to 
determine  the  material  parameter  A  of  the  general  cre^p 
equation  (V-4 ) . 


expressed  as 


tot 


i  L° 

* 

i=n 


where  N=ll  is  the  number  of  divisions.  It  was  found  that 
this  total  average  strain  rate  differs  from  the  creep  rat» 
calculate  only  from  the  highest  temperature  region  by  less 
then  one  percent.  This  shows  that  using  the  six  inches 
uniform  hot  zone  to  be  the  gauge  length  is  a  v^ry  good 
approx imat ion  . 

The  activation  energy  A,‘l  of  the  yarn  was  determined  by 
the  unispecimen  change  of  temperature  technique  which  has 
already  been  discussed  earlier.  The  results  are  shown  in 
Fig. 24  through  29.  The  data  points  in  the  cigutes  were 
fitted  with  a  straight  line  using  the  method  of  linear 
regression.  The  slopes  of  these  lines  are  proportional  to 
the  activation  energy  of  the  impregnated  yarns.  The 
numerical  results  of  the  activation  energy  are  in 
surprisingly  good  agreement.  The  mean  value  for  AH  is  hl.tp 
kcal/mole  with  a  standard  diviation  of  only  2.2  kcal/mole.* 

The  stress  dependence  n  of  the  yarn  w»re  to  h,. 

determined  by  the  unispecimen  change  of  stress  method,  hut 

was  forced  to  change  the  plan  when  unexpected  problems 

arose.  The  first  disadvantage  of  the  uni  specimen  change  o ■ 

stress  method  is  the  changing  of  the  stress,  especially  a: 

*The  physical  interpretations  of  the  activation  energy  ate 
presented  in  Addendum  II. 


VI -C  Creep  Ex per iment :  Results 


The  creep  rate  of  an  impregnated  HM-30D0  fiber  was 
obtained  by  converting  the  rate  of  voltage  change  recorded 
by  the  HP9815  computer  into  elongation  rate  in  mm/hr.  All 
the  temperature  data  points  obtained  during  a  temperature 
raise  were  neglected,  this  is  because  elongation  recorded 
during  this  period  were  resulted  mainly  from  thermal 
expansion  of  the  fiber.  Only  those  datas  points  which  were 
recorded  after  the  furnace  temperature  has  stablized  were 
selected.  Assumption  has  been  made  that  all  of  the  creep 
occurred  at  the  six  inches  uniform  hot  zone  located  at  the 
center  portion  of  the  three  feet  long  furnace.  Analysis 

performed  by  Feldman  [7],  proved  that  this  is  a  very  good 

* 

approx imation . 

In  the  analysis  the  temperature  distribution  of  the 

furnace  was  assumed  to  be  hyperbolic  Fig.  22.  Because  of 

symmetry  of  the  distribution  only  one  side  of  the  furnace 

will  have  to  be  analyized.  The  furnace  was  then 

hypothetically  divided  into  ten  sections  along  its  length, 

and  the  fibers  were  allowed  to  creep  independently  with 

respect  to  the  average  temperature  at  their  own  specific 

location.  A  creep  strain  versus  location  curve  can  then  be 

plotted,  and  the  results  are  shown  in  Fig. 73,  When  the 

average  creep  rate  at  each  section  are  summed  numerically, 

the  result  is  the  total  average  creep  strain  and  can  be 

*  Addendum  I  presents  further  discussion  of  the  gage  length 
and  the  calculation  of  the  pre-exponential  factor  A  of 
Equation  (VI-1) . 
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temperature  in  the  excess  o£  2300  C.  At  this  temperature  any 
change  in  load  can  trigger  failure  even  if  it  was  performed 
by  a  very  steady  hand.  Another  problem  associated  with  the 
changing  of  weiqht  is  that  when  a  load  is  added  or  taken  off 
from  a  fiber,  the  fiber  will  shift,  to  the  right  or  to  the 
left  depends  on  which  action  was  undertaken.  The  sudden 
shift  of  a  fiber  will  interrupt  the  delicately  focused  IRCON 
infrared  pyrometer  which  controls  the  furnace  current.  If  a 
cold  spot  is  exposed  to  the  pyrometer  the  furnace  current 
will  shoot  up  suddenly  and  vise  versa.  It  will  take  a  long 
time  for  the  furnace  to  stabilized  and  even  more  difficult 
to  get  it  back  to  the  exact  temperature  as  it  was  before  the 
load  change.  The  most  curious  phenomenon  was  that  after  a 
fiber  was  hot  stretched  at  a  certain  load  no  creep  was 
detected  at  that  same  temperature  even  when  the  load  was 
increased  by  300%.  It  was  noticed  that  creeping  reassumed 
when  temperature  is  raised  aboved  the  previous  hot 
stretching  temperature,  but  at  a  rate  much  lower  than 
expected  as  shown  in  Fig. 28  and  29.  All  these  factors  made 

it  very  difficult  to  use  the  unispecimen  change  of  stress 
method  to  determine  the  stress  dependence. 

An  alternate  method  was  to  compare  the  creep  rates  at 
the  same  temperature  for  two  constant  load  trails.  One  can 

argue  that  the  unispecimen  method  may  be  is  very  important 
for  working  with  materials  like  highly-oriented 
polycrystalline  graphite,  but  not  necessary  so  for  testing 
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fibers.  Since  all  the  specimens  used  were  made  from  the  same 
roll  of  fiber,  properties  homogeneity  is  aspired. 
Variations  can  only  came  from  the  impregnation  and 
calcination  procedures.  Every  effort  has  been  made  to  ensure 
that  the  specimens'  preparation  procedures  were  as  uniform 
as  possible.  Random  testing  has  confirmed  that  the  matrix 
weight  fractions  were  about  on  various  specimens.  In  the 
light  that  we  have  exactly  identical  dry  fibers,  reasonably 

uniform  impregnation/calcination  procedures  and  also  the 
fact  that  nearly  most  of  the  load  were  carried  by  the  fiber, 
one  can  conclude  that  this  method  is  valid. 

The  calculated  values  for  the  stress  dependence  are 
listed  on  table (3),  and  found  to  be  varied  from  1.2  to  2.1. 
After  careful  eliminations  of  bad  data,  a  mean  n  value  of 
1.41  with  a  standard  deviation  of  .13  was  found.  This 

suggested  that  an  impregnated  HM  3000  fiber  has  some  viscous 
behavior,  for  which  the  creep  rate  is  nearly  proportional  to 
the  applied  stress.  The  stress  dependence  of  viscous  creep 
material  is  one. 
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Oonclus ion 


The  findings  of  the  crack  study  performed  by  Kuhansedgh 
and  Sines  [1]  hinted  that  the  level  of  stress  in  the  fiber 
bundle  may  have  something  to  do  with  the  crack  pattern 
observed.  The  assumption  that  the  radial  fiber  is  the  most 
severely  stressed  fiber  is  confirmed  by  a  simple  analysis. 

The  usefulness  of  a  strong  and  well  bonded  radial  fiber  in 
reducing  processing  anomalies  was  reviewed  by  modifying  the 
computer  model  developed  in  f 1 ] .  Note  that  the  billet  will 

benefit  from  the  radial  fiber  only  if  the  radial  fiber  is 
strong  and  well  bonded.  If  the  radial  fiber  debonds  during 
processing,  its  effectiveness  in  relieving  hoc p  stress  is 

minimal  and  would  be  economically  unfounded  to  weave  in  the 
radial  yarns  to  begin  with.  In  the  light  of  this  problem  two 
different  methods  were  devised,  namely,  the  fiber/matrix 

displacement  mismatch  approach  and  the  stress  relaxation 
approach.  The  aim  of  these  two  methods  were  to  develop  a 
heating  schedule  which  takes  advantage  of  the  creep  action 
of  the  fiber  as  a  stress  relieving  agent,  to  keep  the  radial 
fiber  from  debonding  thereby  preserving  its  intended 
usefulness  to  the  billet. *Two  important  parameters  of  the 

general  creep  equation,  the  activation  energy  and  stress 
dependence  for  an  impregnated  HM  3000  fiber  bundle  were 

determined  by  experiment,  they  took  up  a  value  of  fil.fip 
kcal/mcle  and  1.41  respectively. 

*  Further  analysis  and  discussion  of  the  optimum  time-temp¬ 
erature  path  are  presented  in  Addendum  III. 
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Fig.  Variation  of  Fiber  Volume  Fraction  with  Radi 


ig.3  Process  Flow  Chart  of  the  Billet 
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Fig. 5  Diagram  Showing  A  Characteristic  Area  on  the  Billet's 
Cylindrical  Surface 
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Fig. 7  Billet  With  Equally  Spaced  Fiber  and  Matrix.  (  The  dotted  lines 
repersents  the  shape  that  the  billet  takes  just  after  a  rapid 
heat  up  ) 
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Pulling  Force 


Fig. 8  Forces  Acting  at  the  Fiber  Matrix  interface  due  to  Mismatch  in 
Coefficient  of  expansion,  and  Modulus 
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Resin 


Fig. 10  a.  Diagrammatic  representation  of  deformation  around  a 
discontinuous  fiber  embeded  in  a  matrix  subject  to  a 
tensile  load  parallel  to  the  fiber. 

b.  Stress  distribution  at  the  interface  and  in  the  fiber. 
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Fig. lib  A  typical  stress  vs.  t  plot 
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Fig. 12b  A  typical  stress  vs.  t  plot 
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Fig. 14  The  impregnation  fixture 
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Fig. 16  The  calcination  fixture  loaded  with  specimen  ready  for 
the  tube  furnace 


instrumentation  used  in  the  calcination  process 


Fig. 19  High  temperature  furnace  used  in  the  creep  experiment 


Fig. 19  High  temperature  furnace  used  in  the  creep  experiment 
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Fig. 22  Hyperthetical  temperature  distribution  of  the  creep  furnace 
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STRAIN  RATE /MAX.  STRAIN  RATE 


Fig. 23  Normalized  creep  rate  vs.  location  in  the  creep  furnace 
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Fig. 24  Creep  elongation  rate  vs.  inverse  of  absolute  temperature 
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Fig. 27  creep  elongation  rate  vs.  the  inverse  of  absolute  temperature 


Fig. 28  creep  elongation  rate  vs.  the  inverse  of  absolute  temperature 
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Fig. 29  creep  elongation  rate  vs.  the  inverse  of  absolute  temperature 


W/13000  15V 
5/23/-33 
527  9 

30.’  KCflL'MOLE 


t - 1 —  —  “T-  -  T 


CO 

cn 


I  ,  .  ......... - L.  ...  .  .  I 

il»  O  IP 

.•I 

I  I  I 


c  si  inti  -ayy) 
lltTd  -9N0H3  Ml 

67 


C..  C*‘* 
•I 


Fig. 28  creep  elongation  rate  vs.  the  inverse  of  absolute  temperature 
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PROPERTY  HM-  3000  HM-  10,000 


3 

Density  (  Ib/in  ) 

0.067 

0.066 

Strength  fmpreg.  (ksi) 

350 

355 

ea  (  psi) 

54  x  106 

50  x  106 

ET  (  psi) 

1.5  x  106 

3.4  x  106 

Diameter  (  ym) 

7.2 

7.5 

°A  <  '"F> 

★ 

1.39  x  10"6 

aT  (  /°F) 

★ 

7.5  x  10"6 

Table  1  Conditioned  properties  of  fibers  used  in  this  study  at 
room  temperature  from  their  manufacturer's  data 
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Property 

Cool  tar  pitch 

277CP  -  15V 

Density  (  gr/cc) 

1 .35 

cooking  val ue(  %) 

48.2 

Benzene 

Insolubles{  t) 

12-18 

Quinol ine 

Insolubles( S) 

4-8 

Ash(X) 

0.23 

Table  2  Properties  of  277CP-15V  coal  tar  pitch  matrix  from  their 
manufacturer's  data 
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Append i x 


Calculate  the  Restraining  Hydrostatic  Pressures 


Using  the  original  billet  model,  one  can  solve  for  the 

interfacial  pressure  with  no  radial  effect.  With  these 
informations  one  can  find  the  hoop  stress  in  the  inner  and 
outer  most  fiber, 


c T. 

1 


FP 

[MO,  1) 


_  o  In 

L  (1,?) j  Ri 

tf . 

l 


and 


cr  = 


TP 

[  (58,59) 


-  P 


(59, 60) j  o 


where  ,  RQ  ,  t^  ,  t^  are  the  inside  and  outside  radius 
1  i  o 

of  the  billet  and  fiber  thickness  at  the  inside  and  outside 
locations  respectively. 

The  inside  and  outside  circumferential  strain  are 


4. 

1 


<Ef) 


f 


and 


4 

o 


<E£>* 


The  billet  radius  immediately  after  a  rapid  heat  up  is. 
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Append i x 


Calculate  the  Rest r a  in  ing  Hydrostatic  Pressures 

Using  the  original  billet  model,  one  can  solve  for  the 

interfacial  pressure  with  no  radial  effect.  With  these 
informations  one  can  find  the  hoop  stress  in  the  inner  and 
outer  most  fiber. 
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Tr>  _  P  1  D 
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and 
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fp  -  p  1  r 
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where  ,  R&  ,  t  ^  ,  t^  are  the  inside  and  outside  radius 

ci  o 

of  the  billet  and  fiber  thickness  at  the  inside  and  outside 
locations  respectively. 

The  inside  and  outside  circumferential  strain  are 


< .  = 


1  (Ef>a 


and 


o 


<Ef>a 


The  billet  radius  immediately  after  a  rapid  heat  up  is, 


t£  =  0.035  in 
i 

t£  =0.14  in 
o 

(Ef)a=25xl06  psi  at  5000  F 

<xf  =  1.39xl<rV  F 

/vy  =  4000  F 

A£  =  9.43xl0"4  in2 
.  2 

Ac  =  0.004  in 

A  =  0.007  in2 

c 

o 

p  =0.5 


Addenda  I 

Creep  Parameter  -  The  Linear  Factor  A 

by 

George  Sines  and  Pascal  Got  sis 

One  objective  of  the  study  was  to  determine  the  parameters  in  the  creep 
equation  for  an  impregnated  graphite  fiber  bundle.  The  creep  equation  (VI-1) 
given  on  page  29  is 

~  =  Acnexp(-AH/RT) 

Values  of  n  are  given  in  Table  III  and  discussed  on  page  38.  Values  of  AH 

were  found  from  Figures  24  through  29  and  discussed  on  page  36.  The  value  for 

A  requires  special  consideration. 

In  Figures  24  through  29.  the  ordinate  is  the  natural  logarithm  of  the 
displacement  rate  in  millimeters/hr.  The  specimen  is  of  uniform  cross-section 
without  a  reduction  in  area  in  the  gage  length  so  the  effective  gage  length 
must  be  estimated.  It  has  been  estimated  from  the  temperature  profile  in  the 
furnace  that  the  effective  gage  length  is  1S2  mm  (6  inches).  The  data  point 
for  2700 #C  on  Figure  24  will  be  used  for  an  estimate  of  A  because  it  is  on  the 
best-fit  straight  line  for  the  five  data  points  and  is  a  data  point. 

In  obtaining  a  value  of  A,  uncertainty  exists  or  how  to  calculate  the 

stress.  Should  the  net  fiber  area  be  used  on  the  gross  area  of  the 
impregnated  bundle?  Both  are  meaningful. 
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1 


Using  the  nominal  net  area  of  the  bundle  given  by  the  supplier, 

-4  2 

1.83  x  10  in  ,  and  the  load  of  273  grams  (0.6018  lb)  the  stress  is 
3.18  x  10^  psi.  Using  the  mean  value  of  n  =  1.41,  the  mean  value  of  61.7 
kcal/mol  for  the  activation  energy  T  =  2973°K,  and  a  creep  rate 

-3  -3  _3  _i 

4.76  x  10  hr  ,  we  have  A  =  1.76  x  10  (where  creep  rate  is  in  hr  and 
stress  in  psi).  However,  when  we  use  the  gross  area  of  the  impregnated 
bundle,  which,  admittedly,  was  somewhat  irregular  in  diameter  because  of  the 
pitch,  taking  an  average  diameter  of  0.040  in,  the  stress  is  only  478  psi. 

-3 

Now  we  have  an  A  =  25.5  x  10 

In  our  ongoing  studies  we  will  soon  determine  the  effect  of  excess  pitch 
on  the  creep  rate.  If  it  is  negligible,  then  it  would  be  appropriate  to  use 
the  in  situ  cross-sectional  area  of  the  fiber  bundle. 
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Addenda  II 


Comments  on  Activation  Energy  for  Creep 
by 

George  Sines 

The  creep  tests  were  conducted  on  uni-directional  specimens  in  which  the 
PAN  based  fibers,  "Hercules  Magnamite"  HM  3000,  were  graphitized  by  the 
supplier.  The  specimen  was  impregnated  with  coal  tar  pitch  and  carbonized 
according  to  the  procedure  outlined  on  pages  29-31.  Thus,  the  state  of  the 
bundle  would  be  representative  of  the  first  cycle  of  graphitization  of  the 
billet.  The  measured  activation  energy  of  61.7  teal /mol  for  the  bundle  in 
this  state  may  not  be  representative  of  the  bundle  in  later  cycles  of 
impregnation,  carbonization,  and  graphitization,  in  which  the  matrix  consist 
primarily  of  graphite.  Fully  graphitized  specimens  are  being  prepared  to 
determine  the  effect  of  degree  of  graphitization  of  the  matrix  on  the  creep  of 
the  bundle. 

Values  measured  by  Feldman  [II-I]  on  dry  bundles  of  meso-pitch  fibers  of 
two  kinds  are  66  and  109  kcal/mol.  His  data  for  impregnated  bundles  were  too 
scattered  to  yield  an  activation  energy.  An  even  wider  range  have  been 
measured  for  polycrystalline  carbon  graphite,  60-239  kcal/mol  1 1 1—2  J  . 
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Two  activation  energies  may  be  active  in  the  creep  of  graphites — that 
for  the  carbon-carbon  bond  of  85.3  kcal/mol  and  that  for  lattice  diffusion  of 
163  kcal/mol  [II-2]  .  Measurement  of  activation  energy  may  also  be  distorted 
by  the  large  transient  creep  that  occurs  in  graphite. 

Atom  movement  occurs  in  the  gr aph itiz ation  of  carbon  as  well  as  in 
creep.  Two  peaks  have  been  observed  at  100  and  200  kcal/mole  in  the 
graphitization  of  carbon  [II— 3 ]  .  They  are  thought  to  correspond  to  the  same 
two  mechanisms  postulated  for  creep — carbon  bond  formation  and  lattice 
diffusion. 

Pyrolytic  graphite,  because  its  atomic  structure  is  similar  to  that  of 
graphite  fibers  and  to  that  of  the  matrix  sheaf,  could  have  a  creep  activation 
energy  similar  to  that  of  a  graphitized  bundle.  Values  have  been  measured 
from  100  to  114  kcal/mol  [11-41  . 

The  creep  behavior  and  the  strength  of  multi-cycled,  impregnated, 
graphitized  bundles  in  cylindrical ly  wound  billets  might  be  influenced  by 
hot-working  incurred  by  the  thermal  stresses  during  processing.  Pyrolytic 
graphite  has  had  its  strength  increased  from  an  initial  value  of  26  to  a  final 
value  of  87  ksi  at  2200°C  after  it  had  been  plastically  deformed  at  2750°C 
[11-41  . 

These  uncertainties  complicate  the  development  of  optimum  processing  to 
avoid  failure  of  cy 1 indrical ly  wound  billets  because  of  the  strong  dependence 
of  the  analysis  on  the  activation  energy  for  creep.  We  hope  to  obtain  better 
values  on  bundles  processed  in  a  manner  simulating  the  billet  processing. 
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Addenda  III 

Considerations  for  Optimum  Processing 
by 

George  Sines  and  Pascal  Gotsis 

In  our  study  of  crack  morphology  of  cyl indrical ly  wound,  carbon-carbon 
billets,  it  appeared  that  the  radial  fibers  had  cracked  free  from  the  matrix 

[Text  ref.  1].  There  are  reports  that  in  some  billets,  when  this  debonding 

did  not  occur,  the  radial  fibers  fractured.  In  private  communi eta  ions  Ill-ll  , 

[ 1 1—2 ]  it  has  been  stated  that  the  radial s,  even  though  they  are  debonded,  may 

help  prevent  gross  fracture  of  the  billet.  It  is  thought  the  high  diametral 
expansion  of  the  radial s  puts  presure  on  the  debonded  interface  so  that  the 
radials  assume  load  through  friction,  thereby  reducing  the  tension  in  the 
outer  circumferential  fibers.  This  has  never  been  confirmed.  We  plan  to  test 
this  hypothesis  in  our  proposed  experimental  program. 

To  give  a  qualitative  understanding  of  the  shape  of  the  optimum  time- 
temperature  heating  path  to  avoid  radial  fiber-bond  failure  or  fiber  fracture, 
we  treat  a  simplified  worse  case.  Consider  a  cartesian  lay-up  in  which  the 
fibers  under  consideration  are  a  minority  so  that  their  stress  does  not 
appreciably  retard  thermal  strain  of  the  composite  in  their  longitudinal 
direction.  To  avoid  fracture  of  the  fibers  or  their  debonding,  the  difference 
between  the  thermal  strain  of  the  fibers  and  the  composite  must  be  less  than 
the  critical  elastic  strain  to  cause  failure  plus  the  creep  strain  of  the 
f  ibers. 
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o  .  (t  ) 

(®  -  Of)T(t)  +  /  6  dt 

“  f  r  (o) 

This  equation  is  solved  in  the  text  on  page  23  and  the  symbols  are 
defined  there,  (a  is  the  transverse  coefficient  of  the  composite  in  which 
direction  there  are  very  few  fibers) 

t  =  [<aB  =  ®f )  /A  *cnr.t]  Jexp[AH/RT]dT 

o 

The  integral  does  not  exist  in  closed  form  but  had  to  be  evaluated 

numerically  by  the  IMSL  subroutine  library  of  the  IBM  3033  computer,  the 

initial  heating  to  T  without  creep  which  causes  a  fiber  stress  of  a 

O  Cf  It 

which  has  some  margin  of  safety,  can  follow  any  path.  It  is  at  this 
temperature  the  calculated  path  starts. 

The  material  parameters  used  in  the  calculation  are: 

a  =  5.45  x  10-6  (•C)_1 

II 

af  =  7.72  x  10"7  (•C)-1 
E  =  16.2  x  10**  psi 
n  =  1.41 

A  *  1.81  x  10  (stress  in  psi, creep  rate  in 

hr-1) 

A  was  determined  for  the  best  fit  to  the  data  of  Figure  24.  The  stress 
was  based  on  net  fiber  area. 
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In  Figure  I1I-1  the  path  is  given  for  a  . .  =  130  ksi  and  T  =  1715#K. 

cr  it  o 

To  avoid  failure,  one  hour  must  be  taken  to  raise  the  temperature  the  first 
increment  of  200*K.  After  1.1  hours,  the  temperature  can  be  raised  at  a  very 
rapid  rate,  not  needing  more  than  an  additional  0.02  hr  to  reach  any 
temperature  of  interest  for  graphitization. 

In  Figure  III- 2  for  o  . .  =  100  ksi,  the  initial  temperature  is  1319°E, 

Cf  It 

2 

and  2.2  z  10  hr  is  needed  to  raise  the  temperature  200°K  without  failure. 
However,  only  an  additional  IS  hr  are  needed  to  raise  the  temperature  an 
additional  1500°K. 

In  Figure  III-3  for  a  =  140°K,  the  total  time,  after  reaching  a 

Cf  It 

T  =  1847#K,  is  less  than  0.4  hr  to  reach  3000°K. 
o 

These  calculations  for  cartesian  lay-up  are  very  conservative  when 
compared  to  the  behavior  of  a  cylindrical  lay-up  in  which  the  transverse 
(radial)  strain  is  greatly  constrained  by  the  circumferential  fibers. 

However,  they  do  illustrate  the  extreme  sensitvity  of  the  time-temperature 
path  to  the  failure  strength  of  the  radial  fibers.  They  also  show  the  need  to 
slowly  raise  the  temperature  in  the  initial  period  after  Tq  has  been  reached. 

These  calculations  are  applicable  to  cartesian  lay-up.  The  constraint 

of  the  transverse  fibers  would  have  to  modify  the  thermal  strain 

k  o  . .  (a  -  a,)  in  the  initial  equation, 
cnt  mi 

where  the  constraint  factor  k  would  have  to  be  found  from  micromechanics.  In 
this  simple  illustration  the  material  properties  were  considered  to  be 
constant.  More  accurate  analysis  would  have  them  dependent  on  temperature. 
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III-l  Private  Cononica tion  -  Julius  Jortner,  1982 
III-2  Private  Coaauni cation  -  Steven  Evangel  ides,  1982 
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ADDENDUM  IV 


Stresses  in  the  Radial  Bundles  -  Further  Considerations 

by 

George  Sines 


The  model  used  for  the  first-order  approximation  compared  the  axial 
strains  in  the  circumferential  bundles  at  the  inner  and  outer  diameters  of 
the  billet  to  the  average  axial  strain  of  a  radial  bundle  attached  to  these 
extreme  circumferential  bundles.  (Eq.  1-8,  page  9)  It  was  demonstrated  that 
the  stress  in  the  radial  bundles  exceeds  that  in  the  circumferential 
bundles  when  the  billet  is  heated. 

This  model  neglected  the  local  interaction  of  the  radial  bundles  with 
the  thermal  strain  of  the  composite.  Julius  Jortner  recently  made  a  second- 


order  approximation*  which  gives 


(*°  *  V  4  (ca„  '  Ba>tt 

%  -  V  6 

+  (caQ  -  Ba)AT 


Where  B_ 
G 

r 


B  * 

Ge 


average  axial  mechanical  strain  in  a  radial  bundle. 

axial  tensile  mechanical  strain  in  a  circumferential  bundle 
at  the  outer  radius  Rq.  (It  is  assumed  to  be  approximately 
equal  in  magnitute  to  the  compressive  strain  in  the 
circumferential  bundle  at  the  inner  radius  R^) . 


mechanical  circumferential  strain  in  the  composite  at  Rq. 

coefficient  of  thermal  expansion  of  the  composite  in  the 
circumferential  direction. 


*  To  appear  in  a  forthcoming  report  by  Jortner  Research  to  the  ONR. 
(L.H.  Peebles,  Jr.,  ONR  monitor.) 
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axial  coefficient  of  thermal  expansion  of  the  bundle. 


Notice  that  this  analysis  also  predicts  that  the  strain  (stress)  in 
the  radial  bundles  is  always  higher  than  the  peak  strain  (stress)  in  the 
circumferential  bundles. 

Both  models  are  limited  to  giving  the  average  of  the  axial  mechanical 
strain  in  the  radial  bundles;  it  will  be  seen  that  the  peak  strain  in  the 
radial  bundles  will  be  even  higher:  The  radial  composite  stress  (or  strain) 
which  occurs  upon  heating  is  compressive  and  it  is  zero  at  the  outer  and 
at  the  inner  radii  and  has  a  roughly  parabolic  shape  with  the  peak  near  the 
middle  of  the  radial  thickness  (Reference  1,  Figure  23).  Shear  interaction 
of  the  radial  bundles  with  this  compressive  composite  strain  would  decrease 
the  axial  tensile  strain  in  the  bundles  in  the  middle  and  would  therefore 
increase  it  near  the  outer  and  inner  radii  in  order  to  maintain  the  same 
average  strain.  Thus  the  axial  tensile  stress  in  the  radial  bundles  would 
be  higher  near  the  outer  and  inner  radii  and  higher  than  the  average 
predicted  by  both  of  the  above  analyses.  (Jortner  considers  the  above  in 
a  quantitative  way  in  the  draft  of  his  report  in  preparation.) 

Many  billets  have  a  layered  structure  where  the  number  of  radial 
fibers  is  increased  at  the  outer  radius  in  an  attempt  to  make  the  radial 
bundle  density  more  constant.  Neither  analyses  considers  this  effect  or 
the  coincident  effect  of  termination  of  some  fibers  at  the  layer  interfaces. 
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Rolincik 

P. 

Nilaington 

NA 

AVCO  Systees  Division 

Taverns 

Art 

Nilaington 

KA 

B.F.  GOODRICH  RAD  CENTER 

Stover 

E. 

Brecksville 

OK 

BATTELLE  COLUMBUS  LABORATORIES 

Jelinek 

Frank  J. 

Coluabus 

Ohio 

CALIFORNIA  RESEARCH  A  TECHNOLOGY,  INC. 

Kreyenhagen 

X.  N. 

Chat shot th 

CA 

DEFENSE  NUCLEAR  A6ENCY 

Kohler 

». 

Washington 

DC 

DEFENSE  TECHNICAL  INFORMATION  CENTER 

Alexandria 

VA 

DEPARTMENT  OF  THE  NAVY 

Crone 

J. 

Washington 

DC 

EFFECTS  TECHNOLOGY,  INC. 

Adler 

N. 

Santa  Barbara 

CA 

EFFECTS  TECHNOLOGY,  INC. 

6rahaa 

N. 

Santa  Barbara 

CA 

EXXON  ENTERPRISES,  INC. 

Riggs 

D. 

Fountain  Inn 

SC 

FIBER  MATERIALS,  INC. 

Lander 

L.  L. 

liddeford 

ME 

FIBER  MATERIALS,  INC. 

McAllister 

L. 

Bidder ord 

HE 

6A  Technologies 

Engle 

Glen  D. 

San  Diego 

CA 

GENERAL  ELECTRIC  COMPANY 

Frenke 

Robert 

Cincinnati 

ON 

GENERAL  ELECTRIC  COMPANY 

Hall 

Kenneth  J. 

Philadelphia 

PA 

GENERAL  ELECTRIC  SPACE  DIVISION  -RSO 

Gebhardt 

J. 

Philadelphia 

PA 

HAVES 

Pegg 

R. 

Santa  Fe  Springs 

CA 

HERCULES  CORPORATION 

Christensen 

P. 

Magna 

UT 

ZIPCODE 


94042 

94042 

91702 

95811 

95813 

92663 

90009 

90009 

90009 

90009 

90009 

90009 

90009 

90009 

20332 

♦5433 

45433 

45433 

45433 

45433 

45433 

45433 

45433 

45433 

45433 

87117 

93523 

93523 

02172 

27709 

22314 

22314 

01851 

01887 

01887 

44191 

43201 

91311 

20305 

22314 

20360 

93105 

93105 

29644 

04005 

04005 

92131 

45215 

19101 

19101 

90670 

84044 
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ORGANIZATION 

LASTNAME 

FIRSTNAHE 

CITY 

STATE 

ZIPCODE 

HUGO 

Dyson 

L. 

Gardena 

CA 

90249 

I IT  Research  Institute 

Larsen 

D.  C. 

Chicago 

Illinois 

60616 

JET  PROPULSION  LABORATORY 

Kieeel 

N. 

Pasadena 

CA 

91103 

JORTNER  RESEARCH  1  ENGINEERING,  INC. 

Jortnrr 

Julius 

Costa  Mesa 

CA 

92628 

KAISER 

Fischer 

H. 

San  Leandro 

CA 

94577 

LAURENCE  LIVERMORE  LABORATORIES 

Haiaoni 

A. 

Liver sore 

CA 

94550 

LOCKHEED  MISSILES  l  SPACE  COMPANY,  INC. 

Pineli 

Pat  C. 

Palo  Alto 

CA 

94304 

LOCKHEED  MISSILES  AND  SPACE  COMPANY 

Osaka 

N. 

Sunnyvale 

CA 

94088 

MARTIN  MARIETTA  AEROSPACE 

Koo 

F.  H. 

Orlando 

FL 

32855 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

Uhlaan 

D.  R. 

Caabridge 

HA 

02139 

MATERIALS  SCIENCES  CORPORATION 

Kibler 

J.  J. 

Springhouse 

PA 

19477 

MATERIALS  SCIENCES  CORPORATION 

Rosen 

B.  Halter 

Springhouse 

PA 

19477 

MCDONNELL  DOUGLAS  ASTRONAUTICS  CO. 

Greszczuk 

L.  B. 

Huntington  Beach 

CA 

92647 

MCDONNELL  DOUGLAS  ASTRONAUTICS  CO. 

Penton 

A.  P. 

Huntington  Beach 

CA 

92647 

MCDONNELL  D0U6LAS  RESEARCH  LABORATORY 

Hoi  nan 

H. 

St.  Louis 

NO 

63166 

NASA 

Bersch 

Charles  F. 

Nashington 

DC 

20546 

NASA  Langley 

Houston 

R.  J. 

Langley 

VA 

23665 

NASA  Langley 

Ruaaler 

D. 

Langley 

VA 

23665 

NASA  Langley 

Sawyer 

J.  M. 

Langley 

VA 

23665 

NASA  Marshall  Space  Flight  Center 

Powers 

B. 

Huntsville 

AL 

35812 

NAVAL  RESEARCH  LABORATORY 

»» 

Nashington 

DC 

20375 

NAVAL  SEA  SYSTEMS  COMMAND 

Kinna 

N. 

Nashington 

DC 

20302 

NAVAL  NEAPONS  CENTER 

Jeter 

Eduard 

China  Lake 

CA 

93555 

NSNC/HOL 

Kouhek 

F.  J. 

Silver  Spring 

HD 

20910 

NSMC/MOL 

Rowe 

Charles  R. 

Silver  Spring 

HD 

20910 

OFFICE  OF  NAVAL  RESEARCH 

te 

Pasadena 

CA 

91106 

OFFICE  OF  NAVAL  RESEARCH 

Oincss 

A.  H. 

Arlington 

VA 

22217 

OFFICE  OF  NAVAL  RESEARCH 

Kushner 

A.  S. 

Arlington 

VA 

22217 

OFFICE  OF  NAVAL  RESEARCH 

Peebles 

L.  H. 

Arlington 

VA 

22217 

DSD/DDRAE 

Persh 

J. 

Nashington 

DC 

20301 

PDA  Engineering 

Crose 

J.  6. 

Santa  Ana 

CA 

92705 

PDA  Engineering 

Hack 

T.  E. 

Santa  Ana 

CA 

92705 

PDA  Engineering 

Stanton 

E.  L. 

Santa  Ana 

CA 

92705 

PICATINNY  ARSENAL 

Anazione 

A. 

Dover 

NJ 

07801 

PBATT  t  WHITNEY  AIRCRAFT 

Hiller 

Robert  L. 

Nest  Pal  a  Beach 

FL 

33402 

«*"'TT  t  WHITNEY  AIRCRAFT 

Schaid 

Too 

Nest  Pals  Beach 

FL 

33402 

PURDUE  UNIVERSITY 

Sun 

C.  T. 

Nest  Lafayette 

IN 

47907 

PURDUE  UNIVERSITY 

Taylor 

R.  E. 

Nest  Lafayette 

Indiana 

47906 

RENSSELAER  POLYTECHNIC  INSTITUTE 

Biefendorf 

R.  J. 

Troy 

NY 

12181 

ROCKETDYNE  DIVISION,  ROCKMELl  INTERNTL 

Haukinson 

E.  L. 

Canoga  Park 

CA 

91304 

SAHSO/HNNR 

Bailey 

D. 

Norton  AFB 

CA 

92409 

SAHSO/HNNR 

Brocato 

T. 

Norton  AFB 

CA 

92409 

SANDIA  LABORATORIES 

Northrup 

D. 

Albuquerque 

NH 

87185 

SCIENCE  APPLICATIONS,  INC. 

Eitaan 

D. 

Irvine 

CA 

92715 

SCIENCE  APPLICATIONS,  INC. 

Looais 

Nil lard 

Irvine 

CA 

92715 

SOUTHERN  RESEARCH  INSTITUTE 

Koenig 

John 

Birainghaa 

AL 

35205 

SOUTHERN  RESEARCH  INSTITUTE 

Pears 

C. 

Birainghaa 

AL 

35205 

SOUTHERN  RESEARCH  INSTITUTE 

Starrett 

H.  Stuart 

Birainghaa 

AL 

35205 

SOUTNNEST  RESEARCH  INSTITUTE 

Lankford 

J. 

San  Antonio 

TI 

78284 

SOUTHNEST  RESEARCH  INSTITUTE 

Robinson 

Craig 

San  Antonio 

TI 

78282 

STACKPOLE  FIBERS  COMPANY,  INC. 

Fleaing 

6. 

Loarll 

HA 

01852 

STRATEGIC  SYSTEMS  PROJECT  OFFICE  (PH-11 

Kincaid 

J. 

Nashington 

DC 

20376 
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F1RSTNAHE 

CITY 

STATE 

ZIPCODE 

STRATEGIC  SYSTEHS  PROJECT  OFFICE  (PH-1) 

Neingtr 

S. 

Nashiagton 

DC 

20376 

SUPERTEHP  CENTER 

Leeds 

Don 

Santa  Fe  Springs 

CA 

90670 

SYSTEMS,  SCIENCE  AND  SOFTNARE 

Gurtaan 

6. 

La  Jolla 

CA 

92037 

THIOKOL 

Broaan 

G. 

Brighaa  City 

UT 

84302 

TRW  Systeas 

Kotlensky 

He. 

San  Bernardino 

CA 

92402 

U.S.  ENERGY  RESEARCH  BIV.  ADMINISTRATION 

Littaan 

A. 

Nashington 

DC 

20331 

UNION  CARBIDE  CORPORATION 

Boaaan 

J. 

Cleveland 

OH 

44101 

UNION  CARBIDE  CORPORATION 

Taylor 

A. 

Qak  Ridge 

TN 

37B39 

UNITED  TECHNOLOGIES  RESEARCH  CENTER 

Gallaso 

F. 

Hart Ford 

CT 

06100 

UNITED  TECHN0L06IES-CSD 

Ellis 

Russ 

Sunnyvale 

CA 

9408B 

UNIVERSITY  OF  CALIFORNIA 

BatdorF 

S.  B. 

Los  Angeles 

CA 

90024 

UNIVERSITY  OF  CALIFORNIA 

Sines 

George 

Los  Angeles 

CA 

90024 

UNIVERSITY  OF  NASHlNbfON 

Fischbach 

D. 

Seattle 

NA 

98193 

UNIVERSITY  OF  NYOHING 

Adaas 

D.  F. 

Laraaie 

NY 

B2071 

VIRGINIA  POLYTECHNIC  INSTITUTE 

Hasselean 

D.  P.  H. 

Blacksburg 

VA 

24061 

VIRGINIA  POLYTECHNIC  INSTITUTE 

Jones 

Robert  H. 

Blacksburg 

VA 

24061 

V0U6HT  CORPORATION 

Volk 

H. 

Dallas 

TX 

75211 

WILLIAMS  INTERNATIONAL 

Cruzen 

Scott 

Nailed  Lake 

HI 

48088 

92 


